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17 NMR study of q = spin excitations in a nearly ideal S 

antiferromagnet, Sr 2 Cu0 3 , up to 800 K 
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K.R. Thurber, 1 - 2 !] A.W. Hunt, 1 - 2 T. Imai, 1 - 2 and F.C. Chou 2 
Department of Physics 1 and Center for Materials Science and Engineering, 2 M.I.T., Cambridge, MA 02139 

(February 1, 2008) 

We used 17 O NMR to probe the uniform (wavevector q — 0) electron spin excitations up to 800 K 
in Sr2Cu03 and separate the q — from the q = ±— staggered components. Our results support 
the logarithmic decrease of the uniform spin susceptibility below T ~ 0.015J, where J = 2200 K. 
From measurement of the dynamical spin susceptibility for q = by the spin-lattice relaxation rate 
l/Ti, we demonstrate that the q — mode of spin transport is ballistic at the T — limit, but has 
a diffusion-like contribution at finite temperatures even for T <C J. 
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The one-dimensional Heisenberg spin chain has one 
of the simplest Hamiltonians, H — J^2 t Si ■ Sj+i, yet 
our understanding of its fascinating quantum mechanical 
properties is still developing with recent theoretical (l]-||] 
and experimental [p|-|l2|| studies. A recent breakthrough 
in experimental studies of spin chains is the identification 
of a nearly ideal ID S = h Heisenberg antiferromagnet, 
Sr 2 Cu0 3 , by Motoyama et al. || In this system, S = | 
spins reside at Cu sites, and the superexchange interac- 
tion J is mediated by hybridization with the 2p a orbital 
of the in-chain 0(1) site, see Fig. 1(a). Based on the fit 
of the uniform spin susceptibility x'(q — 0) J is esti- 
mated to be J = 2200 ± 200 K g. 

Sr2Cu03 has proven to be an ideal material for the 
experimental studies of S = h Heisenberg spin chain for 
various reasons. First and foremost, weak inter-chain 
couplings make the temperature of the Neel transition 
(Tjv = 5K) jl3| to a three-dimensional long-range or- 
dered state three orders of magnitude smaller than J, 
Tjf = 0.002J. Thus, the spin excitations of the S = | 
Heisenberg spin chain can be probed at unprecedently 
low scales of temperature and energy. The second ma- 
jor advantage of Sr2Cu03 is that 63 Cu NMR is observ- 
able at the magnetic cation site, because the large J sup- 
presses the nuclear relaxation rates. In a series of pub- 
lications, Takigawa et al. reported detailed 63 Cu NMR 
investigations of the low energy spin excitations [fl0| p~2| . 
They successfully tested the theoretical predictions for 
the q — ±~ staggered mode in the scaling limit in- 
cluding the low temperature logarithmic corrections to 
the staggered dynamical susceptibility Q . The third ma- 
jor advantage of Sr2Cu03, although it has never been 
exploited in the earlier NMR works, is the high local 
symmetry of the crystal structure. The Cu-O-Cu chain 
is strictly straight and the in-chain 0(1) site is located 
in the middle of adjacent Cu sites as shown in Fig. 1(a). 
Therefore the staggered components of the magnetic hy- 
perfine fields from Cu electron spins are canceled out 
at the in-chain 0(1) sites. Accordingly, one can probe 
the low energy spin excitations for the q — long wave- 



length mode (see Fig. 1(c)) separately from the staggered 
q = ±1 mode. Thus 17 NMR study of Sr 2 Cu0 3 pro- 
vides us with a unique opportunity to investigate the 
q = spin excitations down to T ~ 0.002 J in the nearly 



ideal S 



Heisenberg spin chain without interference 



by the staggered mode. Unfortunately, the low oxygen 
diffusion rate in Sr 2 Cu03 severely limits the 17 isotope 
enrichment rate, hence the 17 NMR signal intensity. As 
such, no 17 NMR studies have been reported despite the 
rich information expected for the unexplored q = mode. 

In this Letter, we report the first successful 17 NMR 
investigation of Si'2Cu03 single crystals. We accurately 
measured the temperature dependence of the uniform 
spin susceptibility x'(<7 = 0) by NMR Knight shift at 
the in-chain 0(1) site without suffering from the contri- 
bution by free-spins that limits the accuracy of bulk sus- 
ceptibility measurements and 63 Cu NMR at low temper- 
atures. We found that x'(Q = 0) decreases steeply below 
T ~ 0.015J without the signature of three-dimensional 
short range order approaching the Neel state. Our ob- 
servation supports the presence of a logarithmic decrease 
°f x'(l = 0) a t l° w temperatures, but the quantita- 
tive agreement with existing theoretical models Q is 
not good. Based on 17 NMR spin-lattice relaxation 
rate l/Ti, we also test whether the q = mode of spin 
transport is ballistic or diffusive in a S — \ ID Heisen- 
berg chain, a long standing controversy ^||0. Our data 
strongly support the ballistic nature of spin transport at 
the T = limit, but suggest the presence of diffusion at 
finite temperatures. 

A single crystal of Sr2Cu03 was grown in a traveling- 
solvent floating-zone furnace. 17 isotope was enriched 
into crystals by annealing them at 900-1045 C in 17 02 
gas. NMR measurements were conducted by home made 
NMR spectrometers operated typically at 9 Tesla. The 
hyperfine interaction form factor, F, between 17 nuclear 
spins and Cu electron spins is represented as F\(q) = 
2Ccos(^f) and F 2 (q) = D at the 0(1) and 0(2) sites, re- 
spectively lljj . C and D represent the hyperfine coupling 
constant tensor between the observed 17 nuclear spin 
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and the nearest neighbor Cu electron spins, as schemat- 
ically shown in Fig. 1(a). We determined the hyperfine 
coupling tensor based on the standard K — x plot analy- 
sis as 2C a = 45, 2C b = 95, 2C C = 44, D a = 75, D b = 23, 
D c = 14 ± 10 kOc///B. The superscripts represent the 
crystal axes. One kOe//i£ is a hyperfine interaction of 
gj n h = 4.74 x 1CT 9 cV for 17 0. 

In Fig. 2, we present the temperature dependence of 
the NMR Knight shift K a (l) measured at the in-chain 
0(1) site with 9 Tesla of uniform magnetic field applied 
along the a-axis. We emphasize that one cannot achieve 
such extremely high experimental accuracy at the 0(2) 
and Cu sites in low temperatures, because the enhanced 
staggered spin susceptibility near q = ±ir/a causes both 
homogeneous line broadening as well as inhomogeneous 
line broadening arising from defects |i~2). In contrast 
Fi(q = ± — ) = implies that these line-broadening mech- 
anisms are ineffective at the in-chain 0(1) sites, and the 
linewidth remains one to two orders of magnitude nar- 
rower than at the 0(2) and Cu sites. The 0(1) linewidth 
is only 8 kHz at 30 K, compared to 50 kHz for the 0(2) 
linewidth. The NMR shift K a (l) can be expressed as 

K a (l) = ^— X '(q = Q)+K^ v (l). (1) 

where the first term {K^ pin {l) in Fig. 2) arises from the 
uniform spin susceptibility x'il = 0), while the second 
term Ky V (\) = 0.024% represents the small tempera- 
ture independent contribution of the Van-Vleck term as 
determined by the K — x plot analysis. The temperature 
dependence above 30 K of x'(q = 0) can be fitted very 
well with the high temperature theory |L5p|] . Below 30 K 
(~ 0.015 J), x'il — 0) begins to decrease steeply. We may 
compare the steep decrease of x'(<l = 0) with the theoreti- 
cal prediction by Eggert, Affleck, and Takahashi (lj based 
on conformal field theory. They predicted a logarithmic 
term which has a steep decrease with an infinite slope at 
the zero temperature limit, Jn 2 x'{q = 0) ~ 1 + 2 in{T JT) 
with T a ~ 7.7J, as shown in the inset to Figure 2. At the 
qualitative level, our data support the presence of a steep 
decrease of x'(l = 0). However, our results show that the 
theoretically predicted logarithmic term in x'il = 0) docs 
not account for our experimental results at the quantita- 
tive level for T < 15 K. We note that Motoyama et al. 
H earlier concluded that their bulk susceptibility mea- 
surements matched the theoretical result above ~ 10 K. 
However, we emphasize that our 17 NMR Knight shift 
data do not suffer from the Curie term observed in the 
bulk susceptibility data. As already pointed out by Mo- 
toyama et al. M, the necessity of subtracting the Curie 
term ~ 1/T with a divergently large negative tempera- 
ture coefficient from the bulk susceptibility data leaves 
ambiguities in the temperature dependence of the po- 
tential logarithmic term deduced after the subtraction, 
because the latter has a divergently large positive tem- 
perature coefficient. In contrast, K a (l) directly reflects 



the local spin susceptibility x'{Q — 0) of the neighbor- 
ing Cu electron spins and has no such ambiguity. Thus 
the results in Fig. 2 provide an unambiguous experimental 
test for theories. 

Furthermore, we can rule out the possibility || that 
the steep decrease of x'(l = 0) and/or poor agreement 
with theory is caused by the three-dimensional short- 
range order near the Neel transition at Tjv = 5K based 
on the following consideration. The temperature de- 
pendence of the 17 nuclear spin-lattice relaxation rate 
1/Tf(2) at the 0(2) site is presented in Fig.3(a). 1/Tf 
at 0(1) and 0(2) sites may be expressed as 
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where 7„ is the nuclear gyromagnetic ratio of the ob- 
served nucleus (jn = 5.772 MHz/Tcsla for 17 nuclei), 
and w n (~ 7n-ff) is the resonance frequency. Since the 
form factor at the 0(2) site ^2(9) is independent of 
wavevector q, 1/Ti(2) at the 0(2) sites measures the 
wavevector q integral of x"(<?, Lo n ). This represents the 
strength of the low frequency (u> n ) Cu electron spin fluc- 
tuations. The roughly constant value of 1/Ti(2) re- 
flects the cancellation of the two temperature dependent 
parts of eq.||. 1/Ti oc T which offsets the ~ l/T 1 in- 
crease of the q-integral of the staggered spin susceptibil- 
ity x"(q — i~j w„) down to 10 K, as previously reported 
by Takigawa et al. JlCfl and Thurber et al. 16 based on 
63 Cu NMR. We note that the mild increase of 1/Ti (2) 
down to 10 K is consistent with the logarithmic correc- 
tion intrinsic to ID behavior jl(|. Only at the lowest 
temperature measured (4.2 K) does 1/Ti (2) strongly in- 
crease, indicating the three-dimensional order of the Neel 
state at Tjv = 5 K. We emphasize that the effect of three- 
dimensional short range order is not seen in 1/Ti (2) for 
the region 30 K down to 10 K, where a steep decrease 
of x'(q — 0) is evident in our K a {l) data. This implies 
that the decrease of x'(Q = 0) is a property of the ID 
spin chains, not a 3D ordering effect. On the other hand, 
we cannot rule out the possibility that the discrepancy 
between theory and experiment is caused by the finite 
length of the spin chains caused by defects. The finite 
chain length produces a staggered (q = ± — ) spin density 
oscillation in this temperature regime [p^| which might 
cause unknown effects to x'il = 0). 

Next we turn our attention to the temperature depen- 
dence of the long-wavelength q = mode of the low 
frequency Cu electron spin fluctuations. The temper- 
ature dependence of 1/TfT at the in-chain 0(1) and 
apical 0(2) sites is compared in Fig. 4(a). We found 
that l/T 1 a (l)T may be approximated by an empirical 
form, l/T 1 a (l)T = 0.027 + 4.7 x 10~ 4 T sec" 1 ^ 1 at 
low temperatures. l/T 1 a (l)T shows qualitatively dif- 
ferent behavior from l/T 1 a (2)T, because the hyperfine 
form factor F\(q — ±— ) = filters out the contri- 
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bution of the staggered susceptibility and 1/T"(1)T is 
dominated by q ~ (see Fig. 1(c)). According to the 
theoretical prediction by Sachdev B based on quan- 
tum critical scaling at the low temperature limit, the 
q = contribution to 1/T°(1)T may be written as 
1/T»'«=°T = [(2C b ) 2 + (2C c ) 2 ]< ? 2 7^fc B / 7 r 3 J 2 . The un- 
derlying assumption is that the q — spin excitations 
propagate ballistically without damping at low temper- 
atures rather than diffusive transport. By inserting C b , 
C c , and J into the scaling form, we obtain the theo- 
retical estimate of the contribution by the undamped 
ballistic mode, l/T 1 a ' 9=0 (l)T = 0.029 ± 0.006 sec^R" 1 . 
This is in excellent agreement with our experimental zero 
temperature limit, l/Tf(l)T = 0.027 ± 0.004 sec^K" 1 
without any adjustable parameters. We note that the 
contribution by the q = ±^ branch, l/7f =± " (1)T is 
not negligible at T ^ 0, because the form factor Fi(q) 
will be finite for any q ^ ±— . However by integrating 
the wave- vector dependent scaling form of the staggered 
susceptibility 0] convoluted by the form factor |Fi((7)| 2 , 
we estimate 1/T 1 9_±a (1)T as only of the observed 
rate at 300 K, and the percentage contribution decreases 
slowly with decreasing temperature. Thus the scaling es- 
timation of the sum of the two separate modes of contri- 
butions, 1/T 1 9=0 (1)T+ l/T? =± ° (1)T, while providing a 
very good estimate for T = 0, severely underestimates 
our experimental results at any finite temperature, as 
shown by a dashed line in Fig. 4. Moreover, we found that 
1/Ti(l) depends on frequency as shown in Fig. 3(c). 

Thus, there is an additional contribution to the low en- 
ergy spin susceptibility for q ~ that increases strongly 
with increasing temperature and decreasing frequency. 
This suggests that spin diffusion [fl7| is important, even 
for T <C J ■ We measured the frequency dependence of 
1/Ti(l) between H = 7 Tesla (u n = j n H = 40.4 MHz) 
and H = 14 Tesla (w n = 80.8 MHz) at 77 K and 295 
K. The mild frequency dependence of 1/Ti(l) presented 
in Fig. 3(c) is consistent with the dependence |Ts| ] 

expected for the diffusive contribution, 

-Jj— = [{2C b f + (2C C ) 2 ] ll Y} q = 0) (3) 

where w e = gfi B H/h oc u n . The 63 Cu NMR 1/T X ffij 
also had a small frequency dependent component consis- 
tent with 1 / ^/u^n even though Cu NMR is dominated by 
the q = ±~ modes. l/Ti(2)T for the 0(2) sites does 
not have any frequency dependence within error, indi- 
cating that the dominant q = ±— modes are not diffu- 
sive. To further establish the presence of an unexpect- 
edly large diffusive contribution in the S = h Heisen- 
berg spin chain, we also measured l/T^l) in a related 
one-dimensional spin chain system SrCu02 (see Fig. 1(b) 
and Fig. 3(b)). Since the signal intensity of 17 NMR 
is strong in SrCu02 owing to the higher isotope enrich- 
ment rate, we could measure the frequency dependence of 



1/Ti(l) with higher accuracy. Because of the transferred 
hyperfine coupling D) from the adjacent chain, to 

a good approximation 1/Ti(l) in SrCuC>2 is a superpo- 
sition of the contributions from the q — modes and 
l/Tx{2). This explains why 1/T X (1) in SrCu0 2 asymp- 
totes to 1/Ti(2) at the low temperature limit as shown 
in Fig.3(b). The presence of a large contribution with 
1/ ^/uT^-dependence is clearly seen in Fig. 3(c). 

Theoretically, even whether spin diffusion exists for the 
q = mode of the S = | ID Heisenberg spin chain is 
controversial |Bp|-pJ. Spin diffusion has been measured 
by NMR in S = § spin chains pj and a S = 1 Haldane- 
gap system but to our knowledge not for a S = | 
system. Quantum Monte Carlo results |3|j|] indicate a 
strong increase (faster than T) of 1/Tx(q — 0) with tem- 
perature, but do not determine if the q = peak is truly 
diffusive. Fabricius and McCoy || have calculated that 
the frequency dependence may be ~ w -0 ' 3 rather than 
ui~ - 5 , while a finite size scaling analysis suggests ballistic 
behavior [Q. Our 17 \/T\ results clearly show frequency 
dependence, but are not accurate enough to conclusively 
distinguish the exact exponent. If we assume that the 
extra contribution to 1/T\T is genuinely diffusive with 
frequency dependence oj -0 5 , we can estimate the tem- 
perature dependence of the spin diffusion constant D s , 
which is the only unknown parameter in eq.^[ In such a 
scenario, D s ~ 1/T 2 for T <C J as shown in Fig. 4(b). We 
caution however that this estimate of D s is a lower bound 
on the value, since we are assuming that the additional 
contribution to T% is entirely diffusive. 

To conclude, we have successfully separated the q = 
mode in both static and dynamic spin susceptibility in 
a nearly ideal ID S = | Heisenberg antiferromagnet 
material. We unambiguously demonstrated a steep de- 
crease below T — 0.015J of x'(Q — 0)- Measurements 
of the low frequency q ~ dynamic spin susceptibility 
have a T = limit that agrees with purely ballistic spin 
transport. However, with increasing temperature, the 
dynamic spin susceptibility strongly increases. This re- 
sult establishes the presence of non-ballistic behavior at 
finite temperatures, even for T< J. We suggest the in- 
creased dynamic spin susceptibility is from diffusive con- 
tributions, and estimated a lower bound on the diffusion 
constant D s ~ 1/T 2 . Whether these new results can be 
accounted for by the one-dimensional S = j Heisenberg 
model remains to be seen. 
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FIG. 4. (a)l/Ti°T for in-chain O(l) (•), and 0(2) (o), 
and 1/TfT for 63 Cu (o, divided by a factor of 112) 
in Sr2Cu03. The solid line is the best empirical fit 
to 0(1), 1/T 1 I (1)T = 0.027 + 4.7 x 10~ 4 T sec _1 K _1 . 
The dashed line represents the scaling limit estimate of 
l/T 1 a ' 9=0 (l)T + 1/T 1 °' 9=± «(1)T. (b) Temperature depen- 
dence of the spin diffusion constant D 3 deduced from 
l/T 1 a T(l). This deduced value represents a lower bound 
on the spin diffusion constant. The high temperature limit, 
D s = (J/h)^/2iTS{S + l)/3 = 3.6 x 10 14 sec' 1 , is also shown 
by a dashed line. 



FIG. 1. The fundamental building block of the Cu-0 spin 
chain (Cu •, O o) in (a) Sr2Cu03, and (b) a related zig-zag 
spin chain material SrCu02. Arrows with C, D, and D) 
show the path of transferred hyperfme interactions, (c) The 
spin excitation spectrum of S = \ ID Heisenberg antiferro- 
magnetic spin chain. 



FIG. 2. •: Temperature dependence of the 17 NMR 
Knight shift K a (l) measured at the in-chain 0(1) sites in 
Sr2Cu03. Solid curve represents analytic calculations of uni- 
form spin susceptibility x'(<l — 0) by M. Takahashi et al. for 
J = 2200 K, 2C a = 45 kOe//i B , and K^ v = 0.024%. Inset: 
The same data in an extended scale. The solid curve is the 
theoretically predicted logarithmic form. 



FIG. 3. (a) and (b): 1/T?(2) in Sr 2 Cu0 3 (•), l/"7\ 6 (2) 
in SrCu0 2 (A), and 1/7^(1) in SrCu0 2 (x). (c) Fre- 
quency dependence of l/Tf(l)T in Sr 2 Cu0 3 at 77K (□) 
and 295 K (o), 1/T?(2)T in Sr 2 Cu0 3 at 295 K (•), and 
(l/Ti b (l)T - l/TZ(2)T) in SrCu0 2 at 295 K (X). The 
l/V^ff = = limit is the theoretical estimate of 

1/T 1 9=0 (1)T + 1/T 1 9=± " (1)T. Solid lines are the best linear fit 
with the theoretical constraint at l/\fH = 0. For a contribu- 
tion from spin diffusion, the slope is proportional to l/\/D s . 
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(a) Sr 2 Cu0 3 



(c) Excitation spectrum 
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